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ABSTRACT 


The  pulsed  propulsive  performance  of  a  low-thrust  subliming  -  solid 
ammonium  carbamate  reaction  jet  is  analyzed  and  compared  with  the  results 
of  laboratory  experiments.  Tests  of  the  15  millipound  thruster  were  con¬ 
ducted  at  sea  level  using  a  1. 68  to  1  expansion  ratio  nozzle  and  at  vacuum 
using  a  48  to  1  expansion  ratio  nozzle.  The  transient  processes,  which 
domina  the  short-pulse  or  limit-cycle  mode  of  thruster  operation,  are 
formulated  and  show  good  correlation  with  the  data.  The  apparatus,  pro¬ 
cedures,  and  techniques  required  to  obtain  accurate  test  results  for  a 
low-thrust,  dynamic  mode  of  operation  are  described.  Impulse  bit  size, 
gas  consumption,  and  specific  impulse  are  characterized  in  terms  of 
thruster  geometry,  gas  properties,  and  command  pulse  width  to  provide 
a  systematic  basis  for  design.  Limitations  and  design  criteria  necessary 
for  successful  spacecraft  installation  are  discussed.  A  comparison  is 
made  between  the  performance  of  the  subliming  solid  reaction  jet  system 
and  the  performance  of  other  propellants. 
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NOMENCLATURE 


A  =  area 

a  j  =  Aoa*(2/V)1/2/2Vc 

a2  =-  (Ata  +  /2Vc)[2/(v  +  1  )]  (>  + 1  ,/2  <V  -  1 ) 

a*  =  acoustic  velocity 

B  =  equivalent  orifice-to-nozzle  area  ratio  =  a^/a^ 

=  nozzle  dis  ;harge  coefficient 

Cf  =  thrust  coefficient 

=  thrust  correlation  factor 

Cv  =  exhaust  velocity  coefficient 

D  =  diameter 

F  =  thrust  level 

g  =  proportionality  constant  in  Newton's  second  law 

h  =  heat  of  sublimation 

I  =  specific  impulse 

-  effective  system  specific  impulse 

I  t  =  total  impulse 

Kj  =  empirical  decay  time  integration  factor 

=  choked  nozzle  flow  factor  =  a*  [ 2/(  Y  +  1)  J  ^  +0/2(Y  - 

K 

r 


empirical  rise  time  integration  factor 


Mach  number 


molecular  weight 
Reynolds  number 
pressure 

universal  gas  constant 

temperature 

volume 

specific  volume 
weight  flow  rate 
weight 

pressure  ratio  transformation  =  (1  -  Pc/P 

(1  -  P  /P  )1/Z 
as 

sublimation  coefficient  =  ratio  of  actual  to  maximum  sub¬ 
limation  rates 

ratio  of  specific  heats 

nozzle  expansion  ratio 

effic  iency 

time 

T  /T 
r  r 

density 

working  stress 

dimensionless  time  factor  =  a^0 
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SUBSCRIPTS: 


ambient 

chamber  or  command 

decay 

nozzle  exit 

final 

initial 

nozzle 

orifice 

subliming  solid  propellant 
reservoir  or  rise 
steady  state 
throat 


I.  INTRODUCTION 


There  are  numerous  applications  for  thrust  systems  in  the  10  ^  to  10  *  lb 
range  in  satellite  missions.  Previous  operational  systems  have  utilized  cold 
gas  as  the  propellant.  Cold  gas  systems,  although  fully  developed  and  opera¬ 
tional,  have  weight  and  reliability  disadvantages  associated  with  the  necessity 
of  storing  a  low  density  gas  at  high  pressure. 

In  a  subliming  solid  reaction  jet  propulsion  system,  the  propellant  is 
stored  as  a  solid  in  a  thin  wall  tank.  The  pressure  in  the  tank  is  regulated 
by  controlling  heat  input  to  the  subliming  solid.  The  remainder  of  the  system 
typically  consists  of  a  filter,  propellant  valve,  and  a  nozzle.  Thrust  is 
produced  by  open'ng  the  valve  so  that  the  subliming  vapor  flows  through  the 
nozzle.  Heat  of  sublimation  is  supplied  to  the  solid  to  replenish  the  vapor 
which  is  used.  Since  combustion  and  high  pressures  are  not  involved  in  the 
use  of  the  subliming  solid,  no  reaction  chambers ,  pressure  regulators,  or 
high  pressure  storage  tanks  are  required.  Although  thrusts  ranging  from 
10  ^  to  10  *  lb  may  be  produced  by  subliming  solid  propulsion  systems,  the 
input  heat  required  for  sublimination  is  a  limiting  factor  at  the  higher  thrust 
levels . 

The  pulsed  oropulsion  performance  of  low  thrust  reaction  jets  using 
various  cold  gases  has  been  previously  reported  (Refs.  1-4).  Although 
steady-state  performance  of  subliming  solid  reaction  jets  was  reported  (Ref.  5), 
an  attempt  to  determine  pulsed  performance  was  unsuccessful  (Ref.  6)  (due  to 
corrosion  and  leakage  of  the  sublimation  chamber);  thus  no  information  on 
the  dynamic  performance  of  a  subliming  solid  reaction  jet  is  presently* 
available. 

In  this  document, performance  and  design  relationships  for  a  low  thrust 
subliming  solid  jet  reaction  system  are  analytically  derived  and  are  correlated 
with  the  results  of  sea  level  and  vacuum  experiments  conducted  in  the 


Aerospace  Guidance  and  Control  Laboratory  using  a  15  mlb  thruster  (shown 
schematically  in  Fig.  1).  From  these  relationship*  a  simplified  means  of 
calculating  dynamic  impulse  bit,  gas  consumption,  and  specific  impulse  is 
developed.  Limitations  and  design  criteria  necessary  for  successful  space¬ 
craft  installation  are  discussed. 
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DETAIL  /,  THRUST  STAND 


THRUST  STAND  (DETAIL  I) 
SUBLIMATION  CHAMBER 
QUICK  DISCONNECT 
PRESSURE  GAGE 
ENVIRONMENTAL  OVEN 
THERMOCOUPLE 
SOLENOID  VALVE 
PRESSURE  TRANSDUCER 
LVDT  POSITION  TRANS¬ 
DUCER  (THRUST  PICKUP) 

THRUST  CHAMBER  NOZZLE 
DAMPING  FLUID 

CANTILEVER  BEAM  AND 
PNEUMATIC  FEED 


Fig.  1.  Schematic  of  Experimental  Apparatus 
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II.  ANALYTICAL  APPROACH 


A.  GENERAL  DISCUSSION 

A  good  subliming  solid  propellant  should  produce  low  molecular  weight 
products,  provide  adequate  pressure  for  thrusting  at  reasonable  temperatures, 
and  should  generally  have  a  low  heat  of  sublimation.  Other  investigations 
(Refs.  5-8)  have  found  that  subliming  ammonium  salts  appear  attractive  for 
propulsion  applications.  In  the  present  study  ammonium  carbamate, 
NH^CO^NH^,  was  selected  as  a  typical  subliming  dissociating  solid  for 
analysis  and  testing.  Is  is  assumed  that  the  dissociation  process  proceeds 
fully  as 


NH4COzNH2  +  h  —  2  NH3  +  C02  (1) 

requiring  a  heat  of  sublimation  h  (Ref.  7)  of  877  Btu/lb.  The  molecular 
weight  m  of  the  products  is  26  and  the  ratio  of  specific  heats  y  is  1.31. 

An  experimental  curve  is  presented  of  the  steady-state  vapor  pressure 
P  produced  by  subliming  dissociating  ammonium  carbamate  as  a  function 
of  the  temperature  of  the  gas  in  the  sublimation  chamber  (see  Fig.  2  ). 

Good  agreement  was  obtained  between  the  data  from  Ref.  5  and  the  data 
measured  in  the  Aerospace  Guidance  and  Control  Laboratory.  It  should  be 
noted  that  about  30  minutes  were  required  for  the  vapor  pressure  to  stabilize 
at  each  test  temperature  shown.  This  equilibrium  pressure  time  lag  is  kinet- 
ically  determined  by  bond  breakage  occurring  at  the  surface  of  the  solid 
crystals.  The  sublimation  kinetics  may  be  characterized  by  a  sublimation 
coefficient  a  which  is  defined  (Ref.  5)  as  the  ratio  of  the  actual  rate  of  sublim¬ 
ation  to  the  maximum  equilibrium  sublimation  rate  per  unit  of  solid  surface 

area  at  a  specified  temperature.  It  is  necessary  to  know  both  P  and  a  in 

ps 
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STEADY-STATE  PRESSURE,  Pp,(psio) 


Fig.  Z.  Equilibrium  Vapor  Pressure  Subliming  Ammonium  Carbamate 
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order  to  size  the  total  propellant  surface  area  A^  so  that  the  design  thrust 
level  is  not  kinetically  limited.  ^ 

B.  PRESSURE  TRANSIENTS 


The  reaction  jet  performance  is  based  on  a  dynamical  analysis  of  the 
thrust  chamber  pressure  history.  In  general,  the  chamber  pressure  history 
is  composed  of  steady-state  and  transient  (startup  and  shutdown)  phases 
illustrated  in  Fig.  3.  In  the  following  analysis,  it  is  shown  that  the  transient 
phases  cause  performance  nonlinearities.  Design  relationships  for  pre¬ 
dicting  and  minimizing  these  nonlinearities  are  developed. 

During  startup,  the  flow  rate  through  the  solenoid  supply  valve  orifice 
Wq  exceeds  that  through  the  discharge  nozzle  W^,  causing  the  chamber  pres¬ 
sure  to  buPd  up.  The  expression  reported  (Ref.  4)  for  the  compressible  isen- 
tropic  pressure  buildup  can  be  extended  to  include  the  nonadiabatic  case  where 
the  thrust  chamber  temperature  Tc  does  not  equal  the  reservoir  temperature 
Tf,  If  4  ~  Tc/Tr»  then  the  compressible  expression  for  the  nonadiabatic 
pressure  buildup  is 


dP 

c 

d0 


2a,  t  P 


/  \'n 

/P  /F 

>  \(v+l)/\ 

( v  1 

(is)  -(; 

A 

V  -  7 

V  r/  V 

\) 

nl/2 


t*zt,/zpc  in 


For  steady-state  operation,  dP  /d6  =0,  and  P  /P  =  P  /P  <1.  This  ratio  is 
1  r  c  c  r  s  r 

a  design  variable  and  is  determined  by  the  selection  of  Aq  and  A^. 

Equation  (1)  cannot  be  solved  directly,  but  may  be  simplified  using  one 
of  the  following  assumptions:  (1)  that  the  thruster  nozzle  flow  is  zero  during 
startup;  (2)  that  both  nozzles  are  always  choked;  (3)  that  the  flow  through  the 
orifice  occurs  at  constant  density.  A  comparison  between  the  results  obtained 
using  each  of  these  three  assumptions  and  measured  data  was  made  (Ref.  4), 


2 


It  can  be  shown  (Ref. 


6)  that  A 

P 


Q I  (P 


ps 


F 


1/2 


4 


li 
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Fig.  3.  Typical  Pulse  Pressure  History 


It  was  concluded  that  the  best  correlation  was  obtained  over  a  wider  range 
of  variables  when  assumption  3  was  used.  Thus,  assuming  constant  density 
orifice  flow,  Eq.  (1)  becomes 


Equation  (2)  can  be  integrated  directly,  producing 


(2) 


1/2  i  /  2  *1/2  n  .  n  . 

24  t  =  ln(-y  -  4  B  y  +  1)  + 


1/2, 


1/2 


1/2 


H—B  !„  -Jl  -  4 _ B  -IjB-M) 

/  2  \  ^  /  2  1/2  \1  /2 

(4B2  +4)  (-yZ  -  4 W  By  +  l) 


(3) 


3  ._2 


Assuming  4B  »  4,  Eq.  (3)  is  solved  explicitly  for  P  /P  to  obtaii. 

C  8 


In  the  laboratory  experiments  the  complete  pneumatic  system,  including 
the  sublimation  chamber,  lines,  valve,  and  thrust  stand,  were  all  electrically 
heated  to  the  same  temperature  so  that  4  =  1.  In  Fig.  4  solu  ions  of  Eq.  (4) 
are  presented  for  yQ  =  0.  87  (sea  level)  to  1 . 0  (vacuum) ,  and  %  ^ZB  =  10. 

3 

This  implies  that  AQ/At  >  1  and  affects  the  accuracy  of  Eq.  (3)  in  the  region 
where  P  /P  approaches  unity,  cf.  Refs.  3  anc  4. 

C  6 
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Fig.  4.  Chamber  Pressure  Rise  Transients 


Note  that  the  time  required  to  r«^ih  steady-state  operation  (P^/P^  =  1)  is 

reduced  by  designing  the  jet  reaction  chamber  with  a  small  chamber  volume 

V  and  a  high  orifice  to  nozzle  area  ratio  A  / A,, 
c  0  o  t 

From  Eq.  (4)  the  rise  time  0^,  or  the  time  required  for  the  pressure  in 
the  thrust  chamber  to  reach  steady-state  (P_/Pg  =  !)•  is  expressed  by 


0 

r 


In 


£,/2B 

i rr 


y  +  c, 
7  o  ” 


(5) 


During  shutdown,  the  orifice  valve  is  closed  and  the  gas  accumulated 
in  the  thrust  chamber  discharges  through  the  nozzle.  The  instantaneous 
chamber  pressure  during  an  isentropic  shutdown  was  given  (Ref.  4)  as 

=  -a2(l  -y)0]^/(,^)  (6) 

s 

with  the  cor r mg  decay  time  (time  required  tc  discharge  from  P  to  P  ) 

6  2L 

given  by 

eC=[<Pa'Ps'<,'1')/2''-l]/aZ('  <7> 

However,  for  an  isothermal  expansion  and  decay  (which  might  occur  using  a 
heated  thrust  chamber;. 


dP 

■ar  =  2aa 


(8) 
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JTVI 


so  that  Eq  (6)  becomes 


and  Eq.  (7)  becomes 


p  /p 

c  s 


2a24 


ln(P  /P  ) 
a  s 


A  comparison  between  the  isentropic  and  the  isothermal  chamber  pressure 
decay  transients  is  shown  in  Fig.  5.  It  is  noted  that  there  is  very  little 
difference  between  the  isothermal  and  the  isentropic  decay  pressures.  Since 
one-dimensional  isentropic  flow  is  a  '•o;"'enient  standard  for  comparing 
nozzle  performance,  this  assumption  was  used. 

C.  SYSTEM  CHARACTERISTICS 

The  thrust  and  nozzle  flow  transients  are  readily  determined  from  the 
chamber  pressure  history.  The  dynamic  impulse  bit  size  (total  impulse  per 
pulse)  and  the  dynamic  gas  consumption  (gas  consumed  per  pulse)  are  obtained 
by  integrating  the  rise  transient,  the  steady-state  value  (if  >  0^),  and  the 
d^cay  transient  of  thrust  and  of  gas  ilow.  The  thrust  correlation  factors  and 
nczzle  discharge  coefficients  listed  in  Table  I  (ar.d  discussed  under  losses) 
were  used  in  the  trans;ent  analysis  to  calculate  impulse  bit  size  and  gas  con¬ 
sumption  as  functions  of  command  pulse  0^  (defined  as  the  time  the  valve  is 
open).  The  results  are  shown  in  Figs.  6  and  7,  both  at  sea  level  and  in  a 
vacuum.  Note  that  i.i  the  region  of  small  command  pulses  (0 ^  <  0^),  both 
impulse  and  gas  consumption  are  nonlinear  functions  of  command  pulse.  This 
nonlinearity,  or  deviation  from  the  ideal  square  pulse-wave  (instantaneous 
rise  and  decay),  is  produced  by  the  rise  and  decay  transients.  For  minimum 
imp  ilse  bit  size  (gas  consumpt  on)  at  any  given  command  pulse,  these  transients 
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Fig.  5.  Chamber  Pressure  Decay  Transients 


Fig.  6.  Dynamic  Impulse  Bit  Correlation 


SUBLIMING  AMMONIUM  CARBAMATE 


ynamic  Gas  Consumption  Correlation 


Table  1.  Subliming  Solid  Performance  Characteristics 


Sea  Level  Vacuum 


Average  Gas  Temperature,  *F 

250 

300 

Average  Gas  Pressure,  psia 

60 

45 

Nozzle  Exit  Pressure1 

6  psia 

2480p 

Nozzle  Exit  Temperature4,  *F 

-11 

-310 

Reynolds  Number  at  Nozzle  Throat 

28,  000 

19,  000 

Effective  Nozzle  Expansion  Ratio 

1.6 

45 

Ideal  Specific  Impulse,  sec 

55 

99 

Nozzle  Residence  Time,  psec 

0.  5 

8 

Measured  Steady-State  Thrust,  Fs,  lb 

0.  0145 

0.  0170 

Measured  Steady-State  Weight  Flow,  W8,  10  *  lb/sec 

2.  78 

1.  72 

Thrust  Correlation  Factor,  Ct,  % 

87 

77 

Discharge  Coefficient,  C^,  % 

90 

78 

Impulse  Efficiency,  Cv,  % 

95 

99 

Measured  Specific  Impulse,  sec 

Estimated  Experimental  Accuracy,  % 

52 

96 

@0c  =  20  msec 

5.  0 

@0  =  200  msec 

c 

00 

a 


Isentropic  perfect  gc»8  expansion 
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should  be  minimized  or  possibly  designed  to  effectively  complement  each 
other;  i.  e.  ,  the  sum  of  the  effects  of  the  decay  transient  and  of  the  rise 
transient  approaches  the  ideal  square  pulse-wave.  Equations  (4)  and  (6) 
provide  the  designer  with  relationships  Between  the  variables  so  that  such 
design  tradeoffs  may  be  made. 

Pulsed  specific  impulse  is  shown  in  Fig.  8  as  a  function  of  0c-  For 
vacuum  operation  with  command  mlses  greater  than  20msec,  specific  impulse 
is  independent  of  command  pulse  duration,  although  at  sea  level  a  nonlinearity 
is  produced  by  the  effects  of  transient  operation  at  non-optimum  nozzle 
pressure  ratios. 

D.  PERFORMANCE  LOSSES 

The  tests  in  this  study  were  made  using  a  15  deg  half  angle  conical 
aluminum  nozzle  with  a  measured  throat  diameter  of  0.  188  in.  Two  nozzle 
expansion  ratios  were  used.  The  sea  level  tests  were  performed  using  a 
1. 68  to  1  expansion  ratio  nozzle  and  the  vacuum  tests  were  performed  using 
a  48  to  1  expansion  ratio  nozzle  The  sublimation  chamber  pressure  at  sea 
level  ranged  from  55  to  65  psia  and  was  controlled  by  the  heat  input  from  the 
250SF  environmenta?  oven  and  by  the  discharge  demand  based  on  the  thruster 
duty  cycle.  The  steady-state  vapor  pressures  shown  in  Fig.  2  could  not  be 
maintained  under  flowing  conditions.  As  previously  mentioned,  the  rate  of 
sublimation  is  determined  by  the  amount  of  propellant  surface  area  (provided 
the  sublimation  mechanism  does  not  change  over  the  temperature  range,  e.  g.  , 
surface  spallation).  Since  the  rate  of  vaporization  (or  a)  is  low  (Ref.  5)t  the 
ratio  of  propellant  surface  area  to  nozzle  throat  area  becomes  important  and 

4 

should  be  a  large  number  (typically  >  10  )  in  order  to  minimize  variation  of 
pressure  with  duty  cycle.  In  the  test  setup  the  propellant  surface  area  was 
insufficient  to  maintain  constant  pressure  during  thrusting.  The  vacuum  tests 
were  conducted  using  a  different  duty  cycle  than  at  sea  level  so  that  the  sub¬ 
limation  chamber  pressure  ranged  from  40  to  50  psia  in  a  300°F  environment. 
The  ambient  pressure  in  the  vacuum  chamber  ranged  from  300  to  1000  p 
(depending  also  on  the  thruster  duty  cycle). 
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Perfect  gas  isentropic  flow  relationships  were  used  to  calculate  the 

ideal  performance  characteristics  listed  in  Table  I.  The  throat  Reynolds 

number  (20,  000  to  30,  000)  indicates  that  operation  is  in  the  continuum  flow 

regime  and  that  boundary  layer  viscous  effects  are  not  particularly  severe. 

The  values  of  discharge  coefficient  and  thrust  correlation  factor  are 

based  on  thrust  and  flow  measurements  obtained  in  the  Aerospace  Guidance 

and  Control  Laboratory.  The  sea  level  value  of  C  ,  was  0.  90  and  the  vacuum 

d 

value  of  was  0.  78.  The  sea  level  value  of  was  0.  87  and  the  vacuum  value 
of  Ct  was  0.  77.  For  comparison,  using  unheated  (cold)  gases,  the  average 
value  (Ref.  4)  of  was  0.  85  and  the  average  for  was  0.  80. 

The  isentropic  impulse  efficiency  is  equal  to  the  ratio  of  C^/C^  and 
was  95  and  99  percent,  respectively,  for  sea  level  and  for  vacuum  conditions 

A 

in  the  present  tests.  The  average  value  of  Cv  reported  in  Ref.  4  was  91  per¬ 
cent.  The  high  values  of  C  obtained  during  the  present  tests  using  ammonium 
carbamate  are  attributed  to  the  following  four  effects:  (1)  a  slight  performance 
improvement  (2  to  3  percent)  due  to  some  isothermal  expansion  in  the  heated 
nozzle;  (2)  possible  presence  of  recondensed  ammonium  carbamate  in  the 
pneumatic  system;  (3)  greater  uncertainties  in  the  measurements  due  to 
pressure  and  temperature  fluctuation;  and  (4)  uncertainties  in  the  thermo¬ 
dynamic  properties  of  the  expanding  dissociated  sublimated  gases. 

E.  TRANSIENT  PERFORMANCE  APPROXIMATIONS 

Three  approximate  relationships  were  developed  (Ref.  4)  to  simplify 
preliminary  design  calculations  of  the  dynamic  impulse  bit,  the  dynamic  gas 
consumption,  and  the  pulsed  specific  impulse.  These  expressions  are: 

Impulse  Bit 

For  0  >  0  : 

c  r 

'tot  =  CfA,PslVr  +  (0C  •  V  +  KdV 

(cont.  ) 

3  - 

The  isothermal  impulse  efficiency  was  70%  both  at  sea  level  and  in  vacuum. 


9- 


For  9  <  0  : 
c  r 


‘,0.  ■  CfA.P.<Krer  +  WW 


Gai  Consumption 


For  0  >  0  : 

c  r 


w  =  K  A  P  [K  0  +  (0  -  0  )  +  K,6  .] 

n  t  s  1  r  r  a  r  ddJ 


For  0  <  0  : 

c  r 


w  =  K  A  P  (K  0  +  K  ,0  ,)  (0  /0  ) 

n  t  s  r  r  d  d  r  c 


Specific  Impulse  (vacuum  or  steady-state  at  sea  level) 
I  .  C,/Kn 


where 


K  =  empirical  rise  time  integration  factor 

=  empirical  decay  time  integration  factor 

K  =  choked  nozzle  flow  factor 
n  \ 

=  a*  [2/(Y+  1)](Y  +  1)/2(Y"  U/RT 


It  has  been  found  (Ref.  4)  that  for  vacuum  operation  =  0.  69  and  =  0.  0685. 
The  sea  level  rise  time  integration  factor  is  the  same  as  for  vacuum.  How¬ 
ever,  a  decay  time  integration  factor  =  0.  48  should  be  used  to  provide  good 
correlation  at  sea  level.  The  results  using  this  simplified  technique  of  approx¬ 
imating  the  transient  performance  show  good  correlation  with  the  more  exact 
analysis  and  with  the  measured  performance  (Figs.  6-8). 
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III.  EXPERIMENTAL  PROGRAM 


The  experimental  apparatus  shown  in  Fig.  1  includes  a  gas-generation 
and  pneumatic  system,  a  thrust  stand,  and  an  electronic  command  and  com¬ 
puter  system. 

The  gas -generator  and  pneumatic  system  consists  of  a  sublimation 

3 

chamber  (13.45  in.  capacity),  a  pressure  gauge,  a  thermocouple  temperature 
monitor,  and  a  quick-disconnect  coupler  between  the  sublimation  chamber 
and  the  thrust  stand  assembly.  The  heat  of  sublimation  was  supplied  by 
locating  the  sublimation  chamber  in  an  electrically  heated  environmental 
oven.  Preliminary  tests  indicated  that  the  supply  lines,  solenoid  valve, 
and  thrust  stand  should  be  heated  at  least  to  the  same  temperature  (a  positive 
temperature  gradient)  as  the  sublimation  chamber  for  trouble-free,  non¬ 
clogging  operation.  An  unexplained  hysteresis  between  the  clogging  and 
unclogging  temperature  was  also  observed.  Tests  indicated  that  because  of 
the  limited  vaporization  rate  of  the  solid,  and  differences  in  duty  cycle  and 
in  the  test  configuration,  a  sublimation  chamber  gas  temperature  of  Z50*F 
was  required  at  sea  level  and  of  300° F  in  vacuum  to  establish  the  desired 
operating  pressure  and  to  minimize  fluctuations  during  pulsed  operation. 

A  pressure  regulator  was  not  used  during  the  tests  because  of  materials 
compatibility  and  clogging  problems.  The  pressure  in  the  sublimation  chamber 
was  maintained  at  55  to  65  psia  during  the  sea  level  tests  and  at  40  to  50  psia 
during  the  vacuum  tests  by  proper  adjustment  of  the  time  interval  between 
pulses  (a  Z-sec  interval  for  a  20-msec  command  pulse  and  10  sec  for  ZOO  msec). 
It  is  noted  that  these  sublimation  chamber  pressures  are  lower  than  the 
stabilized  values  shown  in  Fig.  Z  for  the  same  temperature.  The  rate  of  sub¬ 
limation  is  determined  by  the  heat  flux  and  by  the  propellant  surface  area. 

The  long  (0.  5-hr)  stabilization  time  required  to  generate  the  steady-state  sub¬ 
limation  pressures  of  Fig.  Z  was  due  to  insufficient  propellant  area  «n  the 


sublimation  chamber.  Continuous  depletion  of  sublimed  gas  from  the  chamber 
by  the  pulsing  thruster  at  a  rate  greater  than  the  sublimation  rate  therefore 
resulted  in  sublimation  chamber  pressures  which  are  somewhat  lower  than 
the  stabilized  values. 

'  The  complete  apparatus  was  positioned  inside  the  environmental  oven 
during  the  sea  level  tests  so  that  the  lines  and  thruster  stand  were  also  heated 
to  250*F.  In  the  vacuum  tests,  the  gas  flow  lines  and  thrust  stand  were 
heated  by  auxiliary  electric  heaters  to  300*  F,  measured  by  a  bimetallic 
contact -thermometer  mounted  on  the  thrust  stand  in  the  vacuum  chamber. 

The  sublimation  chamber  was  accurately  weighed  before  and  after  a 
test  to  determine  that  ma^s  consumed  during  the  test.  A  tvpical  test  con¬ 
sisted  of  several  thousand  pre-programmed  command  pulse  and  requires 
about  1.  5  hr  to  run.  ^  Each  data  point  is  an  average  of  three  tests. 

The  thrust  measurement  apparatus,  Fig.  9,  has  a  nozzle  mounted  on 
the  end  of  a  short  length  of  tubing  to  serve  the  dual  function  of  a  pneumatic 
feed  line  and  a  cantilever  beam  spring.  When  thrust  is  produced,  the  beam 
is  displaced  an  amount  proportional  to  the  thrust  magnitude.  The  deflection 
of  the  beam  is  sensed  by  a  position  transducer,  and  an  electrical  output 
proportional  to  the  position  (thrust  level)  is  generated.  Viscous  damping 
is  provided  through  the  shearing  action  of  a  flat  plate  moving  through  a 
viscous  fluid.  A  high  natural  frequency  of  the  cantilever  beam  and  nozzle 
assembly  (250  cps)  is  necessary  to  provide  adequate  dynamic  response  for 
the  thrust  transients. 

The  natural  frequency  of  the  thrust  measurement  fixture  was  verified 
experimentally  by  two  methods.  One  method  was  to  displace  the  beam  by  a 
discharge  of  gas  from  the  pneumatic  system  through  the  nozzle.  The  natural 
frequencies  of  the  fixture,  both  with  and  without  damping  fluid,  were  displayed 
on  the  oscilloscope  and  photographed.  The  second  method  was  to  position  a 
coaxial  solenoid  valve  so  that  high  pressure  air  would  discharge  through  the 


An  additional  thermal  stabilization  time  of  about  1.  5  hr  is  required  before 
each  test  to  insure  that  all  components  are  properly  heated 
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valve  and  impinge  on  the  thrust  fixture  nozzle  block.  The  resulting  displace¬ 
ment  of  the  fixture  was  monitored  on  the  oscilloscope.  By  adjusting  the  high 
pressure  air  supply  valve  the  peak  cantilevered  beam  displacement  was  set 
equal  to  the  displacement  experienced  during  thrust  measurement.  The 
thrust  stand  dynamic  response  is  more  than  adequate  for  the  rise  times 
rncountered,  as  evidenced  by  the  correlation  of  thrust  and  pressure^  traces 
in  Fig.  10. 

The  thrust  stand  was  calibrated  in  two  ways.  One  calibration  involved 
the  use  of  a  Scher r-Tumico  force  gauge  and  provided  a  scale  factor  in  terms 
of  volts  per  pound  of  force.  The  other  calibration  was  made  by  orienting  the 
thrust  stand  so  that  the  beam  was  subjected  to  a  1  g  loading.  Test  weights 
were  then  placed  on  the  cantilever  beam,  and  a  correlation  of  force  versus 
electrical  output  (deflection)  was  made.  Both  methods  were  in  excellent 
agreement. 

A  schematic  of  the  command  and  computer  system  is  shown  in  Fip.  11. 
A  preset  counter  or  clock  was  used  to  fix  the  repetition  rate  of  the  command 
pulses,  and  a  digital  counter  recorded  the  number  of  cycles  or  pulses  applied 
to  the  solenoid  valve.  The  output  of  the  thrust  transducer  was  amplified, 
demodulated,  and  fed  to  the  voltage-to-frequency  converter.  The  output  of 
the  converter  was  fed  into  a  reversible  counter  which  is  used  to  integrate 
the  thrust  profile  and  thus  obtain  the  total  impulse  per  pulse  (impulse  bit 
size)  The  reversible  counter  computes  the  impulse  bit  including  drift  o.- 
offset  during  a  time  increment  At.  Then,  the  counter  totals  the  integratec 
offset  and  drift  for  an  equal  time  interval  At  and  obtains  the  difference 
between  the  two  values  to  provide  a  direct  readout  of  total  impulse.  The  use 
of  a  reversible  counter  provides  improved  accuracy  over  the  two  separate 
counters  previously  used  (Ref.  3)  for  the  same  function.  This  is  because 
the  reversible  counter  eliminates  any  errors  due  to  differences  in  timin'’ 
sensitivity  between  individual  counters 


A  pressure  transducer  w 
monitor  pressure. 


a  1200  cps  acoustic  frequency  was  used  io 


The  specific  impulse  was  obtained  by  dividing  the  total  impulse  by 
the  gas  consumed  as  measured  over  a  large  number  of  command  pulses 
(>1000).  The  solenoid  valve  closing  time  was  adjusted  to  equal  the  valve 
opening  time  by  valve  driver  circuitry  modifications  (Ref.  9).  This  con¬ 
figuration  resulted  in  a  valve  open  time  equal  to  the  electrical  command 
pulse  width,  with  a  3.  5  msec  delay  relative  to  the  command  pulse. 

To  complete  the  experimental  evaluation,  the  repeatability  of  impulse 
bit  size  was  examined.  Eighteen  hundred  20-msec  command  pulses  were 
applied  to  the  solenoid  valve,  and  the  total  impulse  of  each  pulse  was  recorded. 
A  attribution  plot  of  these  data  is  shown  in  Fig.  12.  The  spread  and 
character  of  these  data  is  a  result  of  the  thermal  drift  of  the  entire  pneumatic 
assembly.  The  pressure  variations  observed  in  the  supply  tank  during  typical 
operating  conditions  are  due  to  changes  in  sublimation  rate  (caused  by  heat 
flux  variations)  and  directly  affect  the  impulse  bit  size  repeatability. 

An  estimate  of  the  experimental  error  in  measured  specific  impulse 
was  made.  This  estimate  ncluded  the  effects  of:  mass  measurement  errors, 
limited  dynamic  response  of  the  thrust  fixture  ,  calibration  errors,  and 
electronic  errors.  The  largest  error  arises  from  the  uncertainty  in  measur¬ 
ing  propellant  weight.  An  uncertainty  of  2.  5  percent  was  associated  with  the 
weighing  of  4  g  of  propellant  typically  consumed  during  each  specific  impulse 
measurement  consisting  of  2000  command  pulses.  The  second  largest  error 
was  due  to  the  inaccuracy  in  thrust  stand  calibration,  i.  e,  ,  inaccuracies  in 
the  correlation  between  the  displacement  (thrust)  of  the  cantilevered  beam  and 
the  measured  output  voltage.  This  calibration  uncertainty  was  estimated  to 
be  4.  1  percent.  A  third  major  source  of  error  was  produced  by  limitations 
of  the  thrust  stand  dynamics.  This  error  was  estimated  to  be  1.5  percent  at 
a  20  msec  command  pulse  width  and  0.  5  percent  at  a  200  msec  command  pulse 
width.  Electronic  errors  such  as  drift  and  timing  were  insignificant  over 
intervals  of  2000  command  pulses.  A  worst  case  estimate  of  the  average 
experimental  error  was  obtained  by  RSS  (root-sum-square)  of  these  measure¬ 
ment  uncertainties.  As  indicated  in  Table  I,  the  estimated  average  experi¬ 
mental  error  is  5  percent  for  command  pulse  widths  from  20  to  200  msec. 
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IV.  PERFORMANCE  COMPARISONS 


Although  specific  impulse  is  the  most  important  single  measure  of 
propulsive  performance  for  system  design,  consideration  also  must  be 
given  to  the  propellant  density  and  associated  tankage  weight.  The  effective 
system  specific  impulse  is  therefore  defined  (Ref.  10)  as  the  ratio  of  delivered 
total  impulse  to  the  combined  propellant  and  minimum  spherical  tankage 
weight,  or 


where 

P  =  propellant  storage  pressure 

K  =  a  system  factor  to  account  for  electric  heater  and 
minimum  gauge  weights  =  0.  87  for  stored  gases 
and  0.  80  for  stored  liquids  and  solids  (based  on 
actual  spacecraft  designs) 

v  =  the  propellant  specific  volume 

(p/cr)t  =  the  ratio  of  tank  density  to  working  stress  =  3.  25  X  10  ^ft 
for  titanium  with  a  safety  factor  of  2 

Equation  (14)  was  used  to  evaluate  the  effective  system  specific  impulse 
of  H2,  NH^,  N^,  Freon-14,  and  Freon-12  propulsion  systems.  The  conden¬ 
sable  gases  NH^  and  Freon- 12  were  assumed  stored  as  liquids  at  their  room 
temperature  vapor  pressures  of  130  and  85  psia,  respectively.  The  other 
gases  were  assumed  to  be  stored  at  3500  psia.  Solid  ammonium  carbamate 
was  assumed  to  be  heated  by  an  electric  resistance  heater  to  a  temperature 
of  300°F,  producing  a  dissociated  gas  at  a  sublimation  pressure  of  45  psia 
(as  in  the  laboratory  tests).  To  provide  an  equal  basis  for  comparison,  the 
other  five  gases  were  also  assumed  to  be  heated  to  300*F.  No  correction  for 
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possible  partial  isothermal  expansion  effects  {2  to  3  percent  in  impulse)  was 
assumed.  The  results,  listed  in  Table  11  for  a  thrust  level  of  about  10  lb, 
show  that  due  to  its  low  molecular  weight  the  effective  system  specific  impulse 
of  liquid  NH^  is  nearly  20  percent  greater  than  solid  ammonium  carbamate. 
However,  the  system  specific  impulse  of  subliming  ammonium  carbamate 
is  twice  that  of  N^,  Freon- 14,  or  Freon-12  (all  with  I  »  30  sec),  and  5  times 
that  of  (due  to  its  low  bulk  density).  It  should  be  noted  that  while  the  weight 
of  the  heater  was  included  in  this  effective  specific  impulse  comparison,  the 
weight  of  the  electric  pow'  c  source  was  not  included  and  must  be  considered 
in  a  spacecraft  design  application.  It  requires  from  1  to  3  watts /mlb  of  thrust 
to  heat  the  gases  to  300* F.  When  the  heat  of  vaporization  is  included,  about 
10  watts  /mlb  are  needed  for  the  subliming  solid,  and  6  watts  /mlb  for  NH^,  cf. 
Table  II.  If  long  thruster  lines  are  used  downstream  of  the  sublimation 
chamber  (or  boiler),  then  line  heaters  are  required  to  maintain  a  gaseous 
propellant  (to  prevent  flow  clogging),  and  this  additional  heater  and  power 
source  weight  should  also  be  considered.  However,  the  final  selection  of 
one  propellant  over  another  would  depend  on  the  results  of  a  detailed  design 
study  for  a  specific  mission. 
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Table  II  Performance  Comparison  of  Various  Gases 


(T  =  300  *  F,  c  =48,  F  =  15  mlb) 


Gas 

H2 

nh3 

Ratio  of  Specific 

Heats,  y 

1. 40 

1.  31 

Molecular  Weight,  m 

2 

17 

Characteristic 

Velocity,  c*,  fps 

6350 

2220 

Ideal  Specific  Impulse, 
sec 

332 

121 

Delivered  Specific 
Impulse,  seca 

314 

116 

System  Effective 

Specific  Impulse, 
seek 

14 

91 

Power  Required  to 

Heat  to  300° F 
watts /mlbC 

3 

6 

N2 

Freon- 1 4 

F  reon- 1 2 

nh2co2nh 

1. 40 

1.  22 

1.  14 

1.  31 

28 

88 

121 

26 

1690 

1010 

870 

1800 

90 

58 

53 

98 

81 

56 

45 

96 

30 

34 

34 

77 

1 

1.  5 

2 

10 

aMeasured  at  =  70°F  and  adjusted  to  300°F(except  for  NH^CO^NH^) 

^Does  not  include  weight  of  electric  power  supply.  H^,  N£,  and  Freon- 14  are 
assumed  to  be  stored  as  gases  at  3500  psia  with(p^)  =  3.  25  X  10”^  ft  "1  in 
Eq.  (14). 

Q 

Includes  heat  of  vaporization 
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V.  CONCLUSIONS 


Accurate  predictions  of  impulse  bit  size,  gas  consumption,  and  specific 
impulse  can  be  made  using  the  gas  dynamical  relationships  given  for  the 
transient  pressure  histories.  A  simplified  approach  to  transient  performance 
determination  was  also  developed  to  rapidly  evaluate  design  tradeoffs  for 
preliminary  design  purposes. 

Both  the  gas  consumed  and  the  impulse  bit  size  are  nonlinear  functions 
of  command  pulse  width  for  small  pulses  but  become  linear  when  steady-state 
chamber  pressure  is  reached.  The  nonlinearity,  or  deviation  from  the  ideal 
square  pulse-wave,  is  due  to  the  rise  and  dec  y  transients.  For  minimum 
impulse  bit  size  (and  gas  consumption)  at  a  given  command  pulse,  the  rise 
and  decay  transients  should  be  minimized  or  possibly  designed  so  that  the 
sum  of  the  effects  of  the  decay  transient  and  of  the  rise  transient  approaches 
an  ideal  square  pulse-wave. 

The  pulsed  vacuum  specific  impulse  is  independent  of  command  pulse- 
width  for  pulses  >  20  msec  and  is  equal  to  the  steady-state  value.  At  sea  level, 
variable  thrust  coefficients  and  unchoked  nozzle  flows,  caused  by  chamber 
pressure  rise  and  decay  transients,  lower  specific  impulse  and  make  it  a 
nonlinear  function  of  command  pulse. 

No  performance  losses  due  to  the  presence  of  condensed  substances  in 
the  nozzle  flow  could  be  detected.  However,  all  parts  of  the  pneumatic  system 
had  to  be  maintained  at  least  as  hot  as  the  sublimation  chamber,  i.  e.  ,  a 
positive  temperature  gradient,  to  prevent  recondensation  of  the  subliming 
solid  from  clogging  the  system.  A  slight  gain  of  2  to  3  percent  in  specific 
impulse  efficiency  of  the  subliming  solid  compared  to  other  cold  gases  is 
attributed  to  partial  isothermal  expansion  produced  by  this  requirement  for 
a  heated  nozzle. 

The  observed  time  lag  required  to  reach  equilibrium  sublimation 
pressure  is  due  to  kinetic  molecular  reactions  occurring  at  the  solid  surface. 


It  is  necessary  to  know  the  sublimation  rate  (or  coefficient  of  sublimation) 
as  well  as  the  equilibrium  pressure  to  size  the  total  propellant  surface  area 
so  that  sufficient  sublimation  can  occur  to  minimize  large  variations  in 
thrust  level  with  duty  cycle.  This  is  a  significant  factor  to  be  considered 
in  the  design  of  a  subliming  solid  engine. 

The  results  of  a  performance  comparison  between  the  subliming  solid 
system  and  systems  using  five  other  gases  (all  heated  to  the  same  temperature 
and  at  a  thrust  of  10  ^  lb)  show  that  subliming  solids  are  competitive  and 
possibly  even  advantageous  if  sufficient  heating  power  is  available.  The  sub¬ 
liming  solid  heating  power  requirement  of  at  least  10  watts /mlb  of  thrust  is 
mitigated  at  very  low  thrust  levels  or  if  the  duty  cycle  and/or  installation  is 
such  that  a  weight  penalty  for  the  power  source  is  not  incurred.  In  general, 
the  final  selection  of  one  propellant  system  over  another  must  depend  on  the 
results  of  a  detailed  design  study  for  a  specific  application. 
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